Topological flat band models with arbitrary Chern numbers 
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We report the theoretical discovery of a systematic scheme to produce topological flat bands 
(TFBs) with arbitrary Chern numbers. We find that generically a multi-orbital high Chern number 
TFB model can be constructed by considering multi-layer Chern number C — 1 TFB models with 
enhanced translational symmetry. A series of models are presented as examples, including a two- 
band model on a triangular lattice with a Chern number (7 = 3 and an iV-band square lattice model 
with C — N for an arbitrary integer N . In all these models, the flatness ratio for the TFBs is larger 
than 30 and increases with increasing Chern number. In the presence of appropriate inter-particle 
interactions, these models are likely to lead to the formation of novel Abelian and Non-Abelian 
fractional Chern insulators. As a simple example, we test the C — 2 model with hardcore bosons 
at 1/3 filling and an intriguing fractional quantum Hall state is observed. 
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Introduction. The experimental fractional quantum 
Hall effect (FQHE) arises from the highly degenerate 
Landau levels of continuum 2D electron systems, and is 
described by variational wave functions, first proposed by 
Laughlin for the primary FQHE states pQ and later gener- 
alized by Jain for composite fermion states [2] , which are 
analytic functions of the 2D spatial coordinates. Many 
important properties of the FQHE, for example, the hier- 
archy structures and fractionalized excitations [3] H] , can 
be understood within this framework. It even leads to the 
predictions of intriguing non- Abelian fractional quantum 
Hall (FQH) states at certain filling fractions [SHE] . More- 
over, based on a classification of the pattern of zeros of 
symmetric (analytic) polynomials, a systematic way to 
classify FQH states [9j [10] has been proposed. Thus, 
our current theoretical knowledge of FQHE is based on 
the analytic structure the 2D Landau level (LL) Hilbert 
space. 

A key question of fundamental importance in this con- 
text is whether the analytic structure of the candidate 
microscopic wavefunction pQ associated with the LL sub- 
space is essential to the existence of the FQHE phe- 
nomenon. The answer turns out to be no. Conceptually, 
we know that the essential physics of the FQHE can be re- 
garded as the emergence of a macroscopic topological or- 
der [11] at low energy, which is, in principle, independent 
of the microscopic details of the underlying wavefunction. 
Very recently, it has been shown that various Abelian and 
non- Abelian FQH states (also known as fractional Chern 
insulators) can be realized in a large class of so-called 
topological flat band (TFB) lattice models [T2HT4] with- 
out any Landau levels [15-22 ]. The basic idea of the 
TFB models is the following: (a) Use a 2D lattice-based 
tight-binding topological-band model to mimic the non- 
trivial topology of a LL characterized by a nonzero Chern 
number [23], as first proposed by Haldane [24 , and (b) 
find a parameter region to realize a narrow bandwidth 
with a smooth Berry curvature to quench the kinetic en- 
ergy. These new theoretical flatband FQHE discoveries 



not only improve our understanding about the nature 
of the FQHE, but also provide us new ways to realize 
the FQHE in solid state materials [25H29] and ultracold 
atomic systems [I3j [30] . 

Although considerable progress has been made along 
this new direction, the observed FQH states in TFB sys- 
tems have so far been more or less expected since all of 
them can be realized in a LL. This is mainly because the 
Chern number of all these TFB models is C = 1. To ex- 
plore new physics beyond single LLs, it is natural to con- 
sider TFB models with higher Chern numbers [3TJ [32] . 
However, higher Chern number TFBs with large flatness 
ratio (band gap/bandwidth) are much harder to obtain 
without long range hopping. Very recently, a TFB with 
Chern number C = 2 (with a flatness ratio 15) was pro- 
posed [33] . 

In this Rapid Communication, we propose a generic 
and systematic scheme to produce arbitrary Chern num- 
ber TFBs with short-range hopping using multi-orbital 
structures. We believe that these high Chern number 
TFBs will produce novel topological phases when proper 
interactions are introduced. As a simple example, we 
test the C = 2 model with hardcore bosons at 1/3 fill- 
ing and an intriguing FQH state is observed. Moreover, 
TFBs with Chern number C = 2 have the potential to 
realize many new bilayer non-Abelian FQH states [34] 
predicted by the pattern of zeros classification and TFBs 
with Chern numbers C — 8n (n is an integer) may result 
in a class of bosonic integer quantum Hall (IQH) states 
(known as Kitaev's E$ states [35]). 

Generalized Haldane model with Chern number C = 
3. We first consider a triangular lattice model with two 
orbitals per site as shown in Fig. l(a)[ Here, every site 
is colored in red and blue (dark gray and light gray) to 
present the two orbitals A and B. The nearest-neighbor 
(NN) and next-ne xt-ne arest-neigh bor (N NNN) hoppings 
are shown in Fig. l(a)[ while Fig. 1(b) | shows the next- 
nearest-neighbor (NNN) hoppings. The inter-orbital NN 
and NNNN (intra-orbital NNN) hoppings are illustrated 
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FIG. 1: (Color online) The two band model with Chern num- 
ber C = 3 on a triangular lattice, (a) Lattice structure and 
NN (solid bond) and NNNN (dashed curve) hoppings. The 
end point of a link is red (blue) if it is connected to a A (B) or- 
bital, (b) Lattice structure and NNN hoppings. Arrows show 
the sign of the phases. The green disk indicates that each unit 
cell contains one lattice site and the green rectangle represents 
the enlarged unit cell, (c) Single-particle energy spectrum as 
a function of k x and k y . (d) Chiral edge states of the two- 
orbital triangular lattice model, where k± = k x /2 + VSk y /2. 
Each band has Chern number C — ±3. 



by gradient-colored (even-colored) links. The end point 
of a link is red (blue) if it is connected to the A (B) 
orbital. The Hamiltonian of this model is 



(id) {(id)) 
(((id))) 



(1) 



where Al and B\ are the fermion creation operators of 
the two orbitals at site i. The NN, NNN, and NNNN 
bonds are represented by (i, j), ((i, j)), and (((i, j))}. The 
phase factors in the NNN hopping terms are 4>ij — ±tt/2 
with the sign determined by the direction of the ar- 
row [55] . 

In this model, each unit cell conta ins two o rbitals, as 
shown by the green disk in Figs. |l(a)] and |l(b)| and there- 
fore the model has two bands. In order to reach a fiat 
band, we adopt t\ = 1, t<i = 0.39, and £3 = —0.34. The 
NNN hopping is purely imaginary, and hence it is eas- 
ier to realize in optical lattices or solid state materials 
through spin orbit coupl ings. The band structure of this 
model is shown in Fig. l(c)| Due to the particle-hole 
symmetry at (j>ij = ±7r/2, the top and bottom bands 
have the same bandwidth, which is about 1/30 of the 
band gap. Each of the two bands carries Chern number 
±3 [57] • The nontrivial topological structure can also be 
observed by computing the edge states of the system on 
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FIG. 2: (Color online) (a)-(c) Tri- layer Haldane model in an 
ABC-type tri-layer structure. The green rectangle in each 
layer indicates the unit cell of the Haldane model. The dashed 
black circles mark the positions of some sites in the other two 
layers, (d) Chiral edge state of the Haldane model, where 
ki = k x /2 + y/Sky/2. 



a cylinder [Fig. |l(d)] . 

Here we provide an intuitive understanding for this 
C = 3 model and demonstrate its connection to the 
model of Haldane [24 and the ABC-type tri-layer 
graphene structure. We first consider Haldane's model 
on a honeycomb lattice. As shown in Fig. |2(a)| this lat- 
tice has two sites per unit cell, aM and &W, highlighted 
by the green rectangle. The primitive vectors for this 
lattice are = (3^/2,^^/3^0/2) with ao being the 
length of the NN bond. Now we stack three layers of 
this model in a n ABC-type tri-layer pattern as shown in 
Figs. |2(a)||2(c) In the absence of interlayer hopping, all 
three layers share the same Hamiltonian 
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Haldane 



(h3) {((id))) 

t 2 £ c^(aPa™ +6PbW)+H.c > (2) 



where I = 1,2,3 denotes the different layers. It is well- 
known that the lowest band of the Haldane model carries 
a Chern number C = 1 and there exists one chiral edge 
mode at each edge as shown in Fig. |2(d) Therefore, one 



naively expects the tri-layer structure to contain three 
degenerate topological bands with (7 = 1. However, in 
addition to increasing the number of bands, the tri-layer 
structure also enhances the translational symmetry of the 
system. In the absence of inter-layer hopping, the tri- 
layer structure is invariant under the translation by the 
vector (ao, 0) = («i +«2)/3, which is a fraction of the lat- 
tice vector of the original honeycomb lattice, if we also 
permute the layer indices 1 — > 2, 2 — ^ 3, and 3—^1. 
Due to such an enhancement of the translational sym- 
metry, the area of a unit cell for the tri-layer system is 
reduced by a factor of three containing only a single site 
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FIG. 3: (Color online) (a)-(d) Dispersions of the TV-orbital square lattice model along high symmetry directions for TV = 3, 
4, 8, and 20, where t\ — 1, ti — —1/y/N, and = tt/N. In each case the bottom band (red line) is very flat with Chern 
number C — N . (e) The flatness ratio of the bottom flat band increases exponentially with the number of orbit als TV (the 
Chern number TV). 



with two orbitals. As a result, the Brillouin zone is en- 
larged by three times and thus the three-fold degenerate 
(7 = 1 bands form a single band with (7 = 3. It is easy 
to check that this tri-layer model is identical to the tri- 
angular lattice model presented above, if we introduce a 
unitary transformation 



4i 



aW\Al = aM\A]. 



B} 2 



^3 

Bj 3 - 



(3) 



where the subscripts il, i2 and i3 mark different lattice 
sites. 

Generalization to Chern number TV topological flat 
bands. The same idea can also be used to generate 
TFBs with arbitrary Chern numbers. Without loss of 
generality, we consider a square lattice model. We first 
demonstrate this construction using the simplest case 
with Chern number (7 = 2. The lattice and hoppings are 
shown in Figs. Q a) and Qb), where each site contains 
two orbitals (A and B) and the arrows on NN bonds mark 
the directions of positive phase hopping. The Hamilto- 
nian of this model is 



Hl c=2 } =h J2 e ^ A \B 3 + £ t> 3 (A\A 3 + B\B 3 ) 
+ *3 £ (iK+^j+Ec, (4) 

«<M»> 

The NNN hopping amplitudes t\ - are £2 (—£2) along the 
solid (dashed) lines in Fig. [4^b) . Here we adopt the pa- 
rameters £1 = 1, £ 2 = 1/(2 + \/2), £3 = 1/(2 + 2y/2), and 
4>i 3 = ±7r/4 [36 j. In this two-band model, the bottom 
band carries Chern number (7 = 2 and the flatness ra- 
tio is about 30. This Chern number two model can be 
constructed from a bilayer version of the checkerboard 
lattice model first proposed in Ref. [13] (See SI for de- 
tails). Similar to the tri-layer case presented above, this 
bi-layer structure also enhances the translational sym- 
metry and reduces the size of the unit cell by a factor of 
two. 

The generalization to flat band models carrying an ar- 
bitrary Chern number (7 = TV is straightforward. Con- 
sider a square-lattice model with TV orbitals per site. For 



TV > 3, we only need NN and NNN hoppings. The real 
space Hamiltonian is 



^+i]t 



a 



[I] 



i,j 1=1 

+ t 2 [e- < ( 2, - 1 )*Ci3. t lj . +1 +eW- 1 »cM ltj _ 1 



(5) 



The /c-space Hamiltonian is shown in the Supplemental 
Materials. 

For N = 3, the flatness ratio is max imized (~ 58) at 
*2 Ai = -l/VS and (j) = tt/3 [Fig. [5(a) 
case with (7 = N, we adopt t\ — 1, 



For the generic 
£2 -1/VN, and 
= 7r/iV for the A^-orbital model. With these values, 
the flat band is always located at the bottom and the 
Chern number is (7 = N. In the absence of any band 
crossings, one of the higher bands carries Chern number 
C = —(TV - 1)JV, while each of the other (N - 1) bands 
carries Chern number (7 = N. The band structu res o f 
the C — 3, 4, 8, and 20 cases are shown in Figs. |3(a) 



3(d)| We see that the bottom flatband becomes more and 
more flat with an increasing N. As shown in Fig. 3(e)| the 
flatness ratio grows exponentially with the Chern number 
N. The condition £2 = —1/y/N can be further relaxed 
for larger N systems. In particular, it suffices to only 
include NN hoppings when N > 5. 

The TV-orbital (7 = N model can be constructed us- 
ing a TV-layer structure. All the N layers are exactly 
the same up to lattice translations. Each layer is equiva- 
lent to a square lattice subjected to a uniform magnetic 
field. The complex phases of the hopping matrix elements 
are arranged in such a way that each square encircles a 
(j) = 2tt/N flux. Due to the flux structure, the (magnetic) 
unit cell for each layer contains N sites and the lowest 
band of each layer carries Chern number (7=1. If we 
naively put the N layers together, the unit cell has N 2 
sites with the lowest bands TV-fold degenerate. However, 
similar to the models presented above, this TV-layer struc- 
ture has a higher translational symmetry and reduces the 
size of the (magnetic) unit cell by TV times, if we intro- 
duce unitary transformations (simply a permutation of 
TV orbits) similar to Eqs. (|3|. As a result, the Brillouin 
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FIG. 4: (Color online) (a)-(b) The two-orbital Chern number 
two model on a square lattice, (a) NN and NNNN hoppings. 
(b) NNN hoppings. (c)-(e) The 1/3 bosonic FQHE. (c) Low 
energy spectrum E n — E\ versus the momentum k x N y + k y 
for three lattice sizes N s = 18, 24, and 48. (d) Low energy 
spectrum as a function of boundary phase X at a fixed y = 
at 1/3 filling, (e) Spectrum gaps versus 1/N S for three lattice 
sizes. 

zone is enlarged N times and thus the TV-fold degener- 
ate (7 = 1 bands form a single band with C = N. Here 
we emphasize that, although each layer is exposed to a 
magnetic field, the combined Chern number N model is 
realized without an external magnetic field (with zero 
magnetic flux per unit cell and the translational symme- 
try is recovered). 

1/3 bosonic FQH state for the C = 2 model. Similar to 
(7 = 1 TFB models, various FQH states can be stabilized 
in these high Chern number systems once proper interac- 
tions are introduced. As an example, we fill the the (7 = 2 
flatband model [Eq. Q] with hard-core bosons at the fill- 
ing fraction 1/3. For a system with N x x N y unit cells, we 
have N x N y = 3 A/5, where is the number of hard-core 
bosons. The hard-core condition (no more than one par- 
ticle is allowed per site) corresponds to infinite onsite re- 
pulsions: limc/^ UY.i{A A i + B \ B i)( A \ A i + B \ B i- 1 )- 



In Fig.[4](c), we show the low-energy spectra as a func- 
tion of momentum k x N y + k y for different lattices under 
periodic boundary conditions. Here, a sizable spectrum 
gap can be observed, separating the ground state mani- 
fold (GSM) from the other excited states and indicating 
the emergence of an incompressible state. As a definite 
evidence for the emergence of a FQH state, we also com- 
pute the Chern number, which is 2/3 for these three- 
fold nearly-degenerate ground states. The scaling of the 
spectrum gaps is shown in Fig. [4^e) and in Fig. [4] (d), 
we show the low energy spectra as a function of twist 
boundary phase along the x direction with N x = 3 and 
N y = 4. As the phase parameter 6 X changes by 27r, 47r, 
etc., the system goes from one ground state to another, 
and it comes back to the original ground state after X 
varies by 6tt. According to the double layer picture for 
our (7 = 2 model, the most possible candidate for such a 
FQH state is the double layer Halperin (221) state. The 
(7 = 2 model is the simplest model and has the minimal 
flatness ratio in our systematic construction. Therefore, 
we believe FQH states will exist in any of our (7 > 2 
model. 

Conclusion. Before concluding we mention that bilayer 
FQHE is well-known in the 2D LL systems, both exper- 
imentally and theoretically [3&H40] . Multilayer QHE or 
FQHE (beyond bilayers) has not been considered much 
in the continuum systems in contrast to our current work 
where we find that Allayer QHE with a Chern number 
(7 = N should be generic in lattice TFB systems. This 
is potentially an important finding, establishing the TFB 
lattice FQHE to be a much more general theoretical con- 
cept than the ordinary continuum LL-based FQHE so far 
realized in nature in high magnetic fields. 

In conclusion, we propose a generic scheme to con- 
struct TFB models with arbitrary Chern number. These 
new models could be a starting point for many unknown 
Abelian or non-Abelian FQH states with higher Chern 
numbers (> 1) which are generically unknown in the con- 
tinuum 2D FQHE. A more detailed study of these models 
with interactions is beyond the scope of this Rapid Com- 
munication and will be discussed in our future works. 

Note added. After the completion of our work, we no- 
tice that in another recent work [41], a model with high 
Chern number bands was discovered. However, both the 
band gap and the flatness ratio in this model vanish as 
the Chern number increases. After the submission of the 
present Letter, another work [42] observed FQH states in 
our (7 = 3 triangular lattice model, and (7 = 4,5 square 
lattice models. 
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Appendix A: Momentum space Hamiltonian and Berry curvature 

Here we provide some final details for the Chern number C = 3 and C = 2 models proposed in this Letter. We will 
first give out their Hamiltonians in momentum space and then discuss the Berry curvature distribution in momentum 
space. As shown in Figs. 1(a) and 1(b) in the main text, the triangular lattice model is translational invariant along 
the lattice vectors ri = (ao/2, and ?2 — {— &o/2, a/3^0/2), where ao is the length of the NN bond. In 

momentum space, the tight-binding Hamiltonian is 

H [c=s] = ^tn^ c = 3 ^ (Ai) 

k 

where \I/k = (A^B^) is a two component spinor and 7i^ c=3 ^ is a 2 x 2 matrix. Because the intra-orbital (NNN) 
hoppings are purely imaginary, 7i^ c=3 ^ can be expressed in terms of three Pauli matrices, 

^[C=3] = £ M . (A2) 

i 

with i = x, y and z. Here 

h x — — ti [cos ki + cos k 2 + cos(/ci — k 2 )] — ts [cos 2ki + cos 2k 2 + cos(2fci — 2k 2 )] , 

h y = + ti [sin ki — sin k 2 — sin(/ci — k^f\ — ts [sin 2ki — sin 2/^2 — sin(2/ci — 2/C2)] , 

h z = - 2t 2 [sin(fci - 2A: 2 ) - sin(2A:i - k 2 ) + sin(fci + k 2 )} , (A3) 

where k = (k x: k y ), ki = k • ri, and k 2 = k • r 2 . 

In momentum space, the Hamiltonian of the C = 2 square lattice model is given by 

H [c=2] = ^^[ c = 2 ^ k , (A4) 
k 

where = (^k 7 ^k) and 

y[c=2] _ ^ cos 2k x + cos 2/c^) I + V2ti (cos fc^ + cos k y )<7 x — V2ti (cos /c^ — cos k y ) a y — 4t 2 sin k x sin /c^a^. (A5) 

The distributions of the Berry curvature in momentum space for the triangular lattice C = 3 model and the 
square lattice C = 2 model are shown in Figs. g |5 (a) | and S |5(b)| respectively. In both models, the Berry curvature in 
momentum space shows no sharp feature. 



Appendix B: Flatness ratio 

In Fig. |6| we show the contour plot of the ratio of band gap/bandwidth for the nearly flat band at different t 2 and 
ts (20% away from the values we used in the Letter) with ti = 1 and 4>ij — ±tt/2. As can be seen from the figure, in 
a wide region of the parameter space, this ratio remains large indicating a very flat band. 




FIG. 5: (Color online) Distributions of the Berry curvature in momentum space for the flat bands in (a) the triangular lattice 
C — 3 model and (b) the square lattice C — 2 model. 
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FIG. 6: (Color online). The ratio of band gap/bandwidth for the flat band in the triangular lattice C — 3 model. The two 
axes are the NNN and NNNN hopping strengths. The red dot at the center marks the parameters used in the Letter. 




(e) (f) 



FIG. 7: (Color online) The two-orbital Chern number two model on a square lattice, (a) NN and NNNN hoppings. (b) NNN 
hoppings. (c) Band structure as a function of k x and k y . (d) Edge states of this model, (e)-(f) Bi-layer checkerboard lattice 
models. Conventions are the same as in Figs. 1 and 2 of the main text. 



As for the square lattice C = 2 model, its flatness ratio is identical to the checkerboard lattice model proposed 
in Ref. [13]. Therefore, the corresponding contour plot is the same as Fig. l.(b) in the Supplementary Materials of 
Ref. [13]. 



Appendix C: C — 2 square lattice model and bi-layer checkerboard lattice model 



As mentioned in the main text, the two-orbital Chern number = 2 model on a square lattice can be separated into 
two layers of checkerboard models. For the bi-layer checkerboard model shown in Figs. 7(e) and 7(f) t although in 
each layer the unit cell of a checkerboard lattice contains two sites, the unit cell of the bi-layer structure only has one 
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lattice site (with two orbitals). The Hamiltonian of each layer is 

H§ B = tl £ c*«aF 6 J] + £ 4- («F ]t «? + W) +^3 E (<W + « ] ) (CD 

W> «(i,j»> 
where Z = 1,2. By introducing the unitary transformation 

4 = oWt, Aj 2 = a™, 4 = 6^, i4 = 6 [i]t ; (C2 ) 
we recover the two-orbital Chern number = 2 model on a square lattice. 

Appendix D: /c-space Hamiltonian of the iV-orbital C — N model 

The /c-space Hamiltonian of the TV-orbital C = N model reads 

N N 

H [c= N] = j2J2 Mm - mc^cf + [^ m 4' +llt ^ ] + & k + iw 4 +21t 4 + h. c .] , (di) 

k i=i k i=i 

where is the single-particle creation operator of the Ith orbital at momentum k with 1 < I < N and we require 
satisfies periodic boundary conditions, i.e., c\^ +1 ^ = . The functions T and Q are 

Jic(^) = 2t 2 cos(^ + ^ - 
0ik($) = ti[exp(i^) +exp(-zfe 2/ 

^2k(^) = t 2 exp(ife x - iky + z$). (D2) 
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